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Abstract. Resonances or near resonances are ubiquitous among 
the excited nonradial pulsation modes of variable stars and they 
must play an important role in determining their pulsational 
behavior. Here in a first step at nonlinear asteroseismoiogy 
we explore some of the consequences of resonances by means 
of the amplitude equation formalism. 

We show how parity and angular momentum constraints can 
be used to eliminate many of the possible nonlinear resonant 
couplings between modes (and multiplets of modes), and how 
the amplitude equations can thus be simplified. Even when we 
may not be able, nor wish, to make an ah initio computation of 
the values of the coupling coefficients, it is still possible to obtain 
constraints on the nature of the excited modes if a resonance 
between observed frequencies can be identified. 

Resonances can cause nonlinear frequency locking of 
modes. This means that the observed frequencies appear in ex- 
act resonance even though the linear frequencies are only ap- 
proximately in resonance. The nonlinear frequency lock, when 
it occurs, it does so over a range of departures from linear reso- 
nance, and it is accompanied by constant pulsation amplitudes. 
The locked, nonlinear frequencies can differ noticeably from 
their nonresonant counterparts which are usually used in seis- 
mology. This is particularly true for multiplets of modes split 
by rotation. Beyond the regime of the frequency lock, ampli- 
tude and frequency modulations can appear in the pulsations. 
Far from the resonance condition one recovers the regime of 
steady pulsations with nonresonant frequencies for which the 
seismological studies, as they are presently carried out, are justi- 
fied (provided furthermore, of course, that nonlinear frequency 
shifts are negligible). 

Success in identifying a resonance in an observed power 
spectrum depends on the quality of the data. While keeping this 
limitation in mind, we discuss the possible existence of peculiar 
resonances in the pulsations of specific variable white dwarfs 
and 6 Scuti stars. 
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1. Introduction 

Resonances between pulsation frequencies are know-n to play an 
important role in the dynamics of classical variable stars such as, 
for example, Cepheids. (For reviews, see Buchler 1990, 1993.) 
They have a profound effect on the shape of both light and radial 
velocity curves. 

For nonradial pulsators that exhibit a much denser frequency 
spectrum — such as variable white dwarfs, 6 Scuti stars, or y3 
Cepheids — resonances are likely to occur as well. They can 
certainly play a role in modal selection and resonant coupling, 
as has indeed been advocated as the mechanism for saturating 
the amplitudes of 6 Scuti stars (Dziembowski 1993). But reso- 
nances do have other effects that can affect the resulting Fourier 
spectra. Before proceeding, however, we note that resonances 
do not exist between all possible modes in a given star. There- 
fore identifying resonances among observ ed modes may lead to 
constraints on the identification of these modes. 

Resonances couple (nonlinearly) pulsating modes and there- 
fore play an important role in determining the amplitudes of pul- 
sation. The amplitudes of resonating modes are usually smaller 
than would be the case were the modes nonresonantly coupled, 
for instance. In the particular case of a multiplet of modes split 
by rotation, the resonance between the modes of the multiplet 
can be responsible for a strong asymmetric pattern in amplitude 
within a given multiplet (Buchler, Goupil & Serre 1995, here- 
after BGS), as is observed in the case of some variable white 
dwarfs. 

Resonances can also strongly affect mode frequencies and 
the spacings between those frequencies. One way this may hap- 
pen is by “frequency locking” even when linear eigenfrequen- 
cies are not in exact resonance but rather in near-resonance; 
that is, when the exact resonant relation = 0 (kj 

positive or negative integers) between frequencies u>j is not 
exactly satisfied yet nonlinear effects (resonant coupling) can 



160 


J R. Buchler et al.: On the role of resonances in nonradial pulsators 


modify the linear frequencies and cause the observed frequen- 
cies to be locked into exact resonance. As a consequence, a 
near-resonance, through frequency locking, causes the observed 
modes to be equally spaced in frequency. In the case of gravity 
modes (as seen in white dwarfs) this can result in a departure 
from the equal period spacing expected from a first order asymp- 
totic analysis of those modes (see, e.g., Sect. 10.3.3 of Hansen 
and Kawaler 1994). On the other hand, frequency locking within 
a rotationally split multiple! of a rapidly rotating star could yield 
equally spaced frequency splitting, which is to be contrasted to 
the prediction of linear theory where strong departures from 
equal splitting are expected. 

Frequency locking is maintained as long as the modes are 
linearly in near-resonance. This raises the question of how close 
to exact resonance must the linear frequencies be for frequency 
locking to take place. 

It is the purpose of the this paper to address the above issues 
in more detail. 

In order to investigate resonant mode coupling one obvi- 
ously has to go beyond linear effects. Unfortunately, hydrody- 
namical computations are not yet possible for nonradial pul- 
sators. The only theoretical framework presently available is 
the Amplitude Equation (AE) formalism, which has been de- 
scribed in a number of review s (Dziembowski 1993; Buchler 
1988, 1990, 1993) and papers (Dziembowski 1982; Buchler & 
Goupil 1984, hereafter BG84; Goupil & Buchler 1994, here- 
after GB94; van Hoolst 1994). For most well-observed nonra- 
dial pulsators such as white dwarfs, 6 Scuti stars and j3 Cephei 
stars, this theory applies rather well because the theoretically de- 
rived growth rates of the excited modes are very small compared 
to their frequencies. Furthermore, the pulsations are generally 
mildly nonlinear in the sense that the nonlinearities are small 
in the pulsations about the stellar equilibrium configuration. 
For these tw'o reasons the linear and nonlinear periods are very 
close, which by the way is also one of the basic assumptions 
that underlie the current efforts in linear as teroseismology. 

For details of the derivation of amplitude equations in this 
context we refer to GB94 who also give explicit expressions 
for the coefficients of the AEs in an Eulerian formulation. (Van 
Hoolst 1994 has recently given their Lagrangean equivalent.) 
Here we intentionally limit the details of the theory but do in- 
troduce the notation necessary for further discussion. 

The coefficients of the AEs can be computed ab initio, but 
this presents a formidable coding effort which has only been at- 
tempted in a very limited way (Klapp et al. 1985, Dziembowski 
1982), 

It is also possible to extract the values of the coupling coef- 
ficients indirectly from numerical hydrodynamic computations 
of the nonlinear pulsations. This has been done for RR Lyrae 
models in Buchler & Kovdcs (1986a, hereafter BK86a) and for 
Cepheid models in Kovacs & Buchler (1989). It is unfortunate 
that at the present time numerical hydrodynamics calculations 
are not yet feasible for nonradial pulsations in this regard. 

Thus in this paper we take a phenomenological approach 
and undertake a general exploration of the resonant couplings 
that may occur between nonradial modes. 


7 . 7 . Handling of resonances with amplitude equations 

To first order (in amplitude), the velocity pulsations, for exam- 
ple, can be represented by 

u(r,t) = ! ^a u «)e <a -' Q 'L r a (r) (1) 

a 

where a Q and U Q (r) represent, respectively, the complex am- 
plitude and the spatial velocity eigenfunction for the vibrational 
mode a. A similar expression is obtained for the luminosity 
variations, etc. Higher order corrections (BG84, GB94) are 
quadratic, cubic, etc., in the amplitudes {a a (/)} and thus repre- 
sent modal interaction terms. (For completeness we give a brief 
exposition of higher order terms in Sect. 5). The sum is over all 
the dominant modes a (q.v. GB84). We can equivalently count 
a given vibrational mode and its complex conjugate as separate 
modes or add a complex conjugate to expressions such as Eq. 1 . 

If q modes interact to determine the dynamics, then the cor- 
responding AEs have the general form 

^ = ( iu t) +«i)ai +/i(ai,G2, ***,«<?) 

% = (iu> 2 + Kj) a 2 + Zifal, G2 i ' ’ ' > a p) 

# ( 2 ) 

“I" Kq)Oq + fq{0- \ « ^2? ) 

where u; a and K a are, respectively, linear eigenfrequencies and 
growth rates and the functions f a are strictly nonlinear and 
depend on which particular resonances are present. The AEs 
thus form a set of coupled first-order complex nonlinear ordinary 
differential equations, and they govern the temporal behavior of 
the amplitudes and phases. In view of the weak nonlinearity of 
the stellar pulsations we need only consider the lowest relevant 
powers of amplitudes in the functions / Q . In the AEs the 2£ + 1 
components of multiplets (and their complex conjugates) are 
counted as individual modes. Here £ is the angular momentum 
quantum number of the spherical harmonic Yi m (6 , <f>) and 2£ + 1 
corresponds to the number of possible values of m. 

For convenience, one usually converts to real amplitude and 
phase by setting a a = A a e\p(i<p a ). The solutions of the AEs 
come in various forms. On one hand, the (real) amplitudes A a 
can be constant in time, as so-called fixed point (FP) solutions. 
For example, when the AEs have a FP (A a constant and <t> a 
constant) then 

u(rj)= (3) 

Q 

where Cb Q represents the linear frequency corrected for nonlinear 
effects. The pulsation is described to first order by a multi- 
periodic Fourier expansion with q dominant peaks in the Fourier 
spectrum to which should be added “combination frequencies” 
from higher order corrections. 

On the other hand, the amplitudes, A Q , and phases, 0 Q , may 
vary on a long time scale. These modulations can further be 
of several types; namely, periodic, multi-periodic, or irregular 
(chaotic). The phases can, of course, also undergo correspond- 
ing modulations. The periodic or multi-periodic modulations 
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then give rise to peak splittings in the Fourier spectrum when 
the observations span a long enough time interval. 

Many mechanisms can cause amplitude variations. It is well 
known, for example, that the fixed points of 2:1 resonant AEs 
can just be unstable and undergo a Hopf bifurcation leading to 
a periodic variation of the amplitudes and phase. Further away 
from the Hopf bifurcation these modulations can also become 
irregular (see, e.g., Moskalik 1986; BK86a). 

Note, for completeness, that on top of this temporal behav- 
ior one could also have evolutionary changes generally arising 
from, for example, the stellar core’s nuclear evolution (as in 6 
Scuti stars) or from thermal leakage through the star (as in white 
dwarfs). Such secular changes, however, are extremely slow and 
we ignore them here. 

L2. Applicability 

Even though the AEs are easy to write down in general form, we 
caution the reader that a complete description of the pulsations 
in terms of these equations is useful from a practical standpoint 
only when the number of modes is small, or when it is possible 
to consider certain groups of modes in isolation from others. 
Thus, for example, BGS have explored the case of a nonradial 
£ = 1 mode which is rotationally split into three components 
by making the assumption that such multiplets can be studied 
independently of each other. An extension to £ = 2 or higher is 
still straightforward in principle but the number (2£+\) of AEs 
increases rapidly with £ as do the number of resonant coupling 
terms. 

Besides the often large number of excited modes the most 
important puzzle of nonradial stellar pulsations is the frequently 
observed temporal variation of the amplitudes of the excited 
modes. The observations are usually not complete enough how- 
ever to decide whether such variations are periodic, multi- 
periodic, or chaotic, and whether they have a sizeable stochas- 
tic component, or even are stochastically driven. Resonances 
among excited modes can easily cause amplitude modulations. 
In contrast, it turns out to be rather difficult to obtain ampli- 
tude modulations with nonresonant AEs. For example, Buchler 
& Kov^cs (1986b, hereafter BK86b) have shown that in the 
case of only two nonresonantly interacting modes the ampli- 
tudes are necessarily constant. For a larger number of interacting 
modes, modulations do become possible, but stable FPs (with 
fewer modes excited) generally coexist with time-dependent so- 
lutions, and from evolutionary considerations it appears much 
more likely to find the system in one of these FPs. (It would 
take an remarkably strong external kick to make the star jump 
into such a pulsational state.) 

On general grounds we expect the resonant coupling to be 
simpler in white dwarfs than in pulsating Main Sequence stars 
( 6 Scuti and j3 Cephei stars). Indeed, for white dwarfs the pulsa- 
tion frequencies of the #-modes (gravity modes) are very well 
separated from the p-modes (pressure modes) by factors of 1 00. 
Therefore no simple resonances between g- and p-modes oc- 
cur. On the other hand, for the pulsating MS stars the p and 
g spectra have a considerable overlap and many modes have a 


mixed character. The situation here is made even more compli- 
cated by the trapping of p-modes in the interior. Such p-modes 
may thus be active in a resonance, yet their frequencies do not 
appear in the spectrum of the light curve (e.g., Dziembowski & 
Krolikowska 1985). 

Confusing statements have appeared in the literature as to 
the possibility of chaos in stellar pulsators - such as in 6 Scuti 
stars, for example - and we wish to address this point briefly. One 
needs to distinguish between two types of chaos. Rapid, erratic 
amplitude and phase variations, such as those found in R Scuti 
(e.g., Buchler, Serre, Kollath & Mattel 1995, and Buchler et 
al. 1996), which one might term “fast chaos”, require that the 
growth rates of the excited modes be of the same order as the 
frequencies themselves, and thus this regime is clearly beyond 
the range of validity of AEs. On the other hand, for the majority 
of nonradial pulsators (6 Scuti, 3 Cephei, white dwarfs, etc.) the 
growth rates are relatively very small. Chaos, if it occurs in these 
stars, must therefore be in the amplitude and phase modulations 
as a sort of “slow chaos.” A study of such slow chaos is possible 
within the amplitude equation framework (e.g., Wersinger et 
al. 1980; Buchler & Goupil 1988) because k/uj -C 1 is both 
the condition for slow chaos and for the validity of the AE 
formalism. 

In Sect. 2 we examine the general phenomenon of frequency 
locking and its consequences. Interesting, but technical conse- 
quences of frequency locking and its disappearance away from 
the resonance region are given in Sect. 3. In Sect. 4 we con- 
sider resonances where angular frequencies are related by near- 
rational ratios of the form muj a % uuj & or mu; a +71W& +p^ c » 0, 
where m, n and p are (positive or negative) integers. In par- 
ticular, we shall set down selection rules that can be used to 
eliminate various combinations of £ associated with the res- 
onating modes. Sect. 5 is technical and is devoted to proving 
some useful properties of the coupling coefficients. In Sect. 6, 
the possible existence of resonances in nonradial pulsators is 
emphasized with some examples from the literature of variable 
white dwarfs and 6 Scuti stars. With these examples, we then 
illustrate what kinds of constraints could be obtained were a 
resonance to be definitely identified in an observed spectrum. 
We conclude with a summary discussion in Sect. 7. 

The reader who is not interested in technical details is urged 
to skip Sect. 3 and Sect. 5. 

2. Frequency locking - a summary 

A very general result associated with resonances, discussed 
briefly earlier, is the nonlinear phenomenon of frequency lock- 
ing , sometimes also called phase locking . Frequency-lock be- 
tween one two or more resonant modes can occur no matter 
how many additional modes or multiplets may otherwise be 
interacting. 

The reason for this is that for every resonance of the form 
t: kiU)i - o 

i 
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(where the k x are positive or negative integers) we can introduce 
a phase variable of the form = YLki<t>i = 0, and an AE of the 

form d$ a /dt = RHS a , where the RHS a are expressible in terms 
of only the resonant phases (and of course the amplitudes 
A). For the FPs then d$ a /dt = 0 or to the frequency lock 

kiQi = 0 . 

Each resonance can thus give rise to a separate frequency lock. 
2. 7. Two-mode resonances 

We first illustrate frequency locking in the well known case of a 
2: 1 resonance, 2uj a « uf b . Here we consider the situation where 
both modes are linearly unstable, and thus when far from the 
resonance center neither of them becomes a slave mode. (The 
complementary case of a linearly unstable mode a and a linearly 
stable mode b was investigated by Klappet aL 1985, BK86b. In 
that situation the amplitude of the linearly stable mode gradually 
vanishes, more precisely this mode becomes just one of the slave 
modes). 

With the introduction of the real amplitudes and phases, the 
two complex AEs can be rewritten as four coupled AEs in A , 
B , 4> a and <j> b . However, because only the phase combination 
$ = fa - 2$ a appears on the right hand sides, these four AEs 
can further be reduced to the three AEs (cf. Eqs. 18 in Sect. 3) 

dA/dt = K a + F a (A, B, <$>) 

dB/dt =K b + F b (A : B,$) (4) 

d<&/ dt = 6uj + Fj>(A, B , <£). 

The functions F a , F b and F<*> are strictly nonlinear, but they 
are periodic in $ through a cos$ and sin $ dependence. The 
frequency mismatch or “detuning,” which for this resonance is 
SiaJ = u b — 2u; a , appears then in the slow phase. 

When the amplitudes A and B are constant in time, i.e. for 
FPs (dA/dt = dB/dt = 0), the nonlinear frequencies are also 
constant with 

~ u) ■+■ d(p a / dt 
u? b — 2uj + d(j> b / dt . 

Furthermore, d$/dt = d<j> b /dt - 2d4> a /dt=0 and thus $ must 
also be constant. It follows that the observable frequencies obey 

2 u? a = d) b . (6) 

In other words, for stable FP solutions (assuming that they exist) 
the effect of the resonance is then to cause a nonlinear 
synchronization of the two modes. A further consequence 
is that the pulsation exhibits a single period of 27T /£o a (e.g., 
BK86a) so that Eq. (3) reduces to 

u(r, t) = Ae^Uaif) + Be^e i2 ^U b (r) + cc. (7) 

We stress that this synchronization occurs because of the 
nonlinear mode coupling , and despite the fact that the 
linear resonance condition is only approximately satisfied. 


Generalizing, one can now easily verify that a resonance of 
the type nuj a ~ mjj b always introduces a phase 

$ - n<j> a - md b (8) 

and, for FP point solutions (i.e., with A and B constant and 
nonzero), the resonance leads to a frequency lock of the form 

nu) Q = mu) b . (9) 

2.2. Three-mode resonances 

Three-mode resonances are slightly different. A resonance of 
the type muj a & nu? b + p^ c , for example, introduces the phase 

$ = m<p a - nd b — pd> c (10) 

and, when the three amplitudes are constant and nonzero, a 
frequency lock of the form 

mw Q = nuj b + pu; c (11) 

ensues. The general FP solution, instead of involving three in- 
dependent frequencies, can then be expressed in terms of a tw o- 
periodic Fourier expansion with the two basic frequencies h b 
and l «<✓£ » In lowest order 

u(f, t) = Ae i ®e i ^ b+1 ' c)l U a (r) + Be iGjbt V b {r) 

+Ce i ^ t U c (r) + cc. ’ 

Frequency locking can assume yet another form as, for ex- 
ample, in the 1:1:1 resonance within a rotationally-split multi- 
plet considered by BGS. Here, for £ = 1 , it is the phase combi- 
nation 

$ = 0+ + 0_ - 20 o (13) 

that appears in the AEs. The subscripts refer to the three m val- 
ues (-1, 0, +1) of the modes. (We note that this case is mathe- 
matically equivalent to a resonance of the type 2^q % u;+ + .) 

The nonlinear frequency lock, which again occurs when the 
amplitudes are constant, therefore leads to the condition 

u>0 = j(w+ +u)-). (14) 

Or, in other words, the intra-multiplet resonance results in a 
symmetric splitting of the multiplet, even when the linear 
splitting is not symmetric. Note that there is also a nonlinear 
shift in the position of the central frequency peak. 

For an C = 2 mode, an additional phase $ = 0+2 + <£-2 ~ 20o 
appears. For FPs, with all of Aq, Ai~ and A z + constant and 
nonzero, this leads to 

u*> = i(£ +2 + £-2). (15) 

These two peaks are therefore also organized symmetrically 
about the center. 

Thus, quite generally, whatever the resonance may be, one 
always obtains a frequency lock in the vicinity of a reso- 
nance provided a stable resonant FP solution exists . 
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We end this section with a technical note. The previous 
very general result follow s directly from the existence of a cen- 
ter manifold (see, e.g., Guckenheimer & Holmes 1983; Coul- 
let & Spiegel 1984; BG84). Whenever there are one or more 
resonances among the modes {i} of the center manifold — 
that is, relations among the linear frequencies of the form 
53 ki*Ji & 0 , with kj positive or negative integers — then the 
resonant terms in the AEs only involve the specific phase com- 
binations 3> = 53 ki<Pi • The various phase AEs can then always 
be combined into an AE for d$/dt which depends on the am- 
plitudes Ai and $ only. For FPs, when the amplitudes are 
constant, a frequency lock of the form 53 k&i = 0 necessarily 
results. 

23. Condition for frequency locked pulsations 

The above results raise a theoretical question, however. The 
resonant AEs apply not merely near the resonance center, but 
also far from resonance. In contrast, the nonresonant AES hold 
only far from the resonance center. How then are the solutions 
of the resonant AEs related to those of the nonresonant AEs in 
that limit? 

Phrased in terms of observable frequencies, in the limit of 
infinitely large departure from exact resonance for the linear 
eigen-frequencies, one indeed expects to recover a nonresonant 
coupling behavior in w F hich the frequencies are no longer locked. 
What happens to the frequency lock as the linear eigenfrequen- 
cies move away from the resonance and, more specifically, in 
what manner, and how far from the resonance center does fre- 
quency locking break down? 

The question is also of practical interest. First, in view of 
the possible constraints it can bring up (see next section), one 
should wonder whether, in stellar conditions, frequency locks 
can ever occur, and if so, over what range of parameters. A 
second point concerns the special case of a multiplet of modes 
split by rotation. With the rotation considered a perturbation, the 
linear eigenfrequencies of a l = 1 triplet of modes, for example, 
are obtained as usual and schematically as: 

UJ m ~ ^0 + AiJm / 1 s-\ 

= ^od + m£ E, + m 2 (^) 2 H 2 + . . •) 

with Q the rotational frequency and E i , S 2 factors of order unity 
depending on the stellar structure and on the nature of the modes. 
It is then clear that the eigenfrequencies obey the near-resonant 
relation 2u?q = u;+ + + 8u) y with a resonance frequency mis- 
match ~ O(fi/uo) 2 . Relations (16) are commonly used 

to infer rotation rates from the observed splittings. This can he 
dangerous because nonlinear effects can vitiate this kind 
of approach. Indeed, when frequency locking occurs and non- 
linear effects are large enough they can cause equal splitting far 
beyond the small Q regime. However, as we will show here, for 
large enough H, the frequency lock is eventually broken. In this 
regime one expects modulated amplitude and frequency pulsa- 
tions, and as one moves very far from the resonance condition 
the amplitude modulations become increasingly shallow 7 and 


negligible. Here one recovers the regime of pulsations with non- 
resonant, nonlocked frequencies, i.e. with those predicted by the 
nonresonant AEs w'hich differ from the linear frequencies only 
by nonlinear corrections, i.e. terms of order of the amplitudes 
squared. Only in this regime are the current scisniologi- 
cal studies justified as they compare observed frequencies 
with nonresonant ones. (This statement of course further sup- 
poses that the mentioned nonlinear frequency shifts are small). 

We now turn to a more detailed discussion of the solutions of 
the AEs as one moves aw 7 ay from the resonance condition. For 
the situation to be of interest the parameters must be such that ( I ) 
a stable multimode resonant FP exists at the resonance center, 
and (2) the nonresonant AEs (valid ’infinitely' far from reso- 
nance) also have a stable multimode FP (with the same modes). 
This implies in particular that we assume that the relevant modes 
are linearly unstable (Otherwise, far from the resonance, the res- 
onant linearly stable modes join the multitude of slave modes, 
i.e. their amplitudes become negligibly small, as for example in 
the bump Cepheids, cf. BK86a). 

A first important result is that, quite generally, a resonant 
FP solution with frequency locking can exist in the vicinity of 
a resonance (here we summarize the rather technical Sect. 3); 
more specifically, for 

\6u>\/uj&0(k/u;) ( 17 ) 

where k is the growth rate of the mode. This is a rough order-of- 
magnitude estimate with factors of 1 0 possible on either side. As 
is obvious, the nature of the solutions of the AEs depend on the 
numerical values of the linear and nonlinear coefficients. Note 
that rel. (17) holds also for splittings, and whether the cause of 
the splitting is rotational or magnetic. 

In turn, does the absence of frequency locking imply that 
the amplitudes are time-dependent ? The AEs are too com- 
plicated to study analytically, but numerical solutions provide 
some answers as to the nature of the solutions (see Sect. 3). 
One can distinguish roughly three regimes when one moves 
away from the resonance center i.e when Suo/h is increased: 

(1) It often happens that near the resonance center {8^/k small) 
one has a stable FP for which the quoted frequency lock there- 
fore exists. There, the frequency may substantially differ from 
their nonresonant frequency counterparts which are used in seis- 
mological studies. 

(2) As 6jc/k is increased this FP eventually becomes unstable, 
or it disappears. In either case, in this second regime one wit- 
nesses a bifurcation to another solution which we find to be 
of an oscillatory nature in the numerical examples considered, 
thus giving rise to amplitude modulated pulsations. However, as 
Slc/k gets large the amplitude of this modulation decreases and 
eventually evanesces as 1 /8u>. Concomitantly the period of the 
modulation decreases as In/bu). The satellite-peaks (due to the 
amplitude modulation) thus gradually move and their amplitude 
vanishes. 

(3) Finally, one thus approaches the third regime, viz. a ’clean’ 
steady spectrum in which the frequencies are given by their 
(unlocked) nonresonant values. For further details we refer the 
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reader to Sect. 3, which also shows a few typical Fourier spectra 
that one may expect. 

To summarize, near a resonance one has a frequency locked 
spectrum, whereas further away from the resonance one obtains 
peaks at the } (which are no longer locked). These peaks may 
be sharp, or they may be broadenend because of amplitude and 
frequency modulation. 

As an example, consider the effects of rotation. The 6 Scuti 
stars, for instance, typically have large rotational velocities v ~ 
150-200 km/s and stellar radii about twice the solar radius. 
This leads to a large and asymmetric rotational splitting. This 
asymmetry thus gives values of the off-resonance parameter 
6u)/ujo ~ 10" 3 -10~ 2 in the range of observed frequencies. The 
growth rates for these stars are small, k/uj ~ 10' 7 up to 10~ 4 
for the most unstable modes in the observed frequency range 
(Dziembowski 1993). This gives a range of Suj/k of 10 - 10 5 
for S Scuti stars. 

On the other hand, variable white dwarfs are believed to be 
slow rotators and S^j/ujo ~ 10~ 8 - 10~ 6 . Theoretical computa- 
tions of growth rates for DAV stars yield k/uj ~ 10~ 10 — 10 -5 
(Dolez & Vauclair 1981; Lee & Bradley 1 993; Bradley & Winget 
1994). Here the range of 6uj/k is 10“ 3 - 10 +4 . 

Give or take an order of magnitude, the condition for a fre- 
quency lock (Eq. 17) is satisfied and it can occur for the most 
unstable modes and for modes that predominantly sample the 
regions of slower rotation. Since the approximate condition for 
frequency locking depends on linear growth rates which can 
vary strongly from one part of the power spectrum to another, 
frequency locking may well occur only for different subsets of 
modes in a particular star. 

Although growth rates remain very uncertain, the above or- 
ders of magnitude indicate that one can expect frequency lock- 
ing in variable white dwarfs and slowly rotating 6 Scuti stars, 
whereas it is less likely to occur in rapidly rotating <5 Scuti stars. 

3. The extent of frequency locking - theoretical 
considerations 

The main results of this section were already summarized in 
Sect. 2. A reader not interested in theoretical fine-points of fre- 
quency locking may therefore skip this section without preju- 
dice. Here we examine what happens to the frequency locked 
solutions of the AEs as one moves away from the resonance 
center, and we show that these fixed points must eventually dis- 
appear, giving rise to amplitude modulations. 

3.1. 2:1 resonance 

We first consider the case of the 2:1 resonance, to be specific. 
The appropriate resonant AEs are given in Eq. (18) 

^ = K a A - cia a A* - Qa b AB 2 + R a ABcos($ - S a ) (18a) 

at 

^ = K b B - q* b A-B - q*B 2 + R b A 2 cos{$ + 6 b ) (1 8b) 

at 

^7 =6u>- ( q[ b A 2 - q bb B 2 ) + Hq^A 2 + q{ b B 2 ) 
at 


A 2 

+R b sin($ + <5fc)— — 2R a sin(<F - S a )A (18c) 

t> 

where the superscripts R and I refer to the real and imaginary 
parts, respectively, and where r k = R k exp(i6 k ), for k=a,b. 

Consider the third equation (18c). Since the last two terms 
are bounde d ( [s in [ < 1, A and B finite by virtue of Eqs. (18a, b) 
A ~ 0(y/K/q) ^ 0(k/R), one notes that as Suj gets large 
there can no longer be a DM FP solution after some critical 
value 6u> b i /. For the resonant terms to play a role the quadratic 
and cubic terms must contribute to the same order, hence we 
have k m qA 2 ^ RA. It is therefore easy to see that one has 
to have 6uj bl f/uj = 0{k/uj). This is of course an important 
practical result. However, we need to stress that this is an order 
of magnitude result since we have assumed here that the real 
and imaginary parts of the q and r are of the same order. 

We wish to add a technical comment here. Prima facie, the 
conclusions of the previous paragraph appear at odds with those 
of Buchler & Kovacs (1986a) who, in the context of the bump 
Cepheid resonance, found the DM FP solution to exist for a 11 
Suj — ► —oo, with B — * 0. The reason is that different assump- 
tions were made. In Buchler & Kovacs it was assumed that 
K b < 0 so that in the NR limit only a SM FP exists. Far from 
the resonance mode b becomes just another slave mode. Here 
both k > 0, and both amplitudes stay of the same order as Sjj 
is increased. 

When we are far from resonance then the nonresonant (NR) 
AEs should apply. For the case of tw o modes they are given by 

^ - (q* A 2 + q n b B 2 )A ( 1 9a) 

at 

^ = K b B - (q*A 2 + q«B 2 )B (19b) 

at 

Strictly speaking the standard cubic coefficients q jk in 
Eqs. (18) and q jk in (19) differ by terms which are 0(\/6uj) 
(BK86a, Eqs. 12), and thus vanish in the limit of large Suj. The 
question again arises as to how the solutions of the resonant and 
nonresonant AEs compare in the limit of large Suj). 

The resonant AEs cannot be solved analytically. For the 
purpose of illustration we consider a numerical example. We 
stress that the results we obtain depend on the values of the 
coefficients, both quantitatively, but also qualitatively. Thus for 
example the nature of the bifurcations can be affected. 

In order to put ourselves into the situation where a double- 
mode FP exists in the nonresonant limit ( i.e. A and B different 
from zero) we have to assume that both n a and K b are positive, 
and that the discriminant D = gaaQbb ~ QabQba > 0 8(rt0* 

(The case where the a mode is linearly stable and the b mode is 
unstable has been considered in BK86a, and the reverse case in 
Moskalik 8i Buchler 1990). 

In this numerical example we choose the parameters values 
(« a = 1.0, = 2.0, C = 1 .0, q* b = 0.5, = 0.5, q£ b = 4.0, 

R a = 15.0, R b = 3.0, 8 a = 6 b = 0, and all qfj = 0. Note that 
with this choice of 6 phases all qf k = qf k ). For convenience and 
without loss of generality we have effectively rescaled the time 
in the AES by setting k q - 1 . 
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Fig. 1 . 2:1 Resonance: Amplitudes and phase as a function of the 
resonance detuning parameter Su; in units of k q ; solid line : stable and 
dashed: unstable frequency locked solution; dotted line: nonresonant 
FR 


The FPs of the NR AEs and their stability can easily be found 
( cf. Buchler& Kovdcs 1 986). There are 2 unstable SM FPs with 
respective amplitudes, A = 1.0, B = 0.0 and B = 0.7071 , A = 
0.0. Further they have a stable DM FP with amplitudes A = 
0.8944 and B = 0.6325. 

In Fig. 1 we present a bifurcation diagram for the FPs of 
the resonant AEs Eqs. 1 8. Because of our choice of parameters 
the bifurcation diagram is symmetric around 8w = 0. From 
the top down we show the amplitudes A and B and the phase 
4>. The solid lines denote the stable FPs, and the broken lines 
correspond to other unstable FPs. One observes that a stable 
FP point exists only for a range |<5u;| < ~ 18.24. In the 

domain 12.64 < 6u> < 18.24 four FPs are seen to coexist, but 3 
of them are unstable. One of the unstable branches (long dashes) 
exists for all 8uj. Far from the resonance this unstable solution 
approaches the unstable single-mode NR FP A = 1,5 = 0. 
For « 21.16 the A amplitude of one of the unstable FPs 
(short dashes) vanishes, and interestingly at this point the system 
bifurcates to the unstable NR FP A = 0, B = 0.707. 

In the region 8u> > 8 ^ / where no stable FPs exist, a nu- 
merical integration of the AEs leads to an oscillatory solution 
(limit cycle). In the star this corresponds then to periodically 
modulated pulsations. 

A sample for several values of 8u is shown in Fig. 2. In these 
examples the A amplitude is greater than B. The oscillations are 
seen to occur about the double-mode FP of the NR AEs, viz. 
A - 0.8944, B = 0.6325 (shown as thin horizontal lines). For 
large <5u> the the amplitude modulations become sinusoidal and 
their amplitudes are evanescent, scaling with \/8u). The modu- 


lation frequency, on the other hand, becomes very large and ap- 
proaches 8uj. As <5oTis decreased the amplitude modulations first 
reach a plateau at 6u> % 15.5, and then in a very short interval 
their period gets large and they become very distorted. We have 
not found any oscillatory solutions below 8u ; ~ 15.495. (How- 
ever, they can be hard to find and a more thorough search might 
uncover them). Note that in the region 1 5.495 & 8u & 1 8.25 the 
oscillatory solution coexists with the stable phase-locked FP. 

This coexistence of two possible types of pulsational states 
also raises the question of how the star selects its pulsational 
state. This clearly depends on its past history, and on fluctuations 
that are present and can make the star jump from one branch to 
another. 

In this example we have illustrated the correspondence be- 
tween the solutions of the resonant and the NR AEs. The two 
AEs agree in the limit 8uj — > oc where they should, but, short 
of the limit, the solution of the resonant AEs undergoes ex- 
tremely rapid oscillations with vanishing amplitude. Such small 
differences are expected in the reduction to normal forms which 
extract the essence of the pulsations and discard finer details. 

3.2 . Rotationally- split multiplet- 1:1:1 resonance 

We now turn to the 1:1:1 resonance which was discussed in 
connection with a rotationally induced splitting in an V = 1 
multiplet (BGS). Because of the cancellation of many of the 
resonant terms (due to angular momentum constraints) only the 
combination phase $ of Eq. 1 3 appears. Rather than being those 
of two 1:1 resonances, mathematically the AEs are the same 
as one would obtain for a three-mode resonance of the type 
2 u>i^ lu '2 + u? 3 . This is the reason why only the parameter 8uj = 
oj + + tu— — 2 u ; 0 rather than the separate frequency differences 

— cuo and u^o — appear. 

For a triplet £ = 1 of modes split by rotation of respective fre- 
quencies , o; 0 , the corresponding amplitudes A- , A+, Aq 
and phases <£ 0 , 0+ then obey the resonant AEs (BGS): 



= k_j4_ + R-AlA + cos ($-6_) 
- A-(q u A 2 _+ q n A\ + q^Al) 

( 20 a) 

Ao 

— kqAq + RqAqA^A— cos(4>+<5o) 



— Ao(q2)A 2 _ + <722^5 + to-4;) 

( 20 b) 

A + 

= k+A + + R+A\A„ cos ($-8+) 
- A+(qu AL + < 732 + < 733 - 4 +) 

( 20 c) 

$ 

= 8u> - 2R$A- A+ sin(4>+<5o) 



sin($~<5_) + sin(4>-<5+)j 

( 20 d) 


The quantity 8u = u;+ + a>_ — 2u;o measures the departure from 
exact resonance which, here, is also the departure from equal 
frequency splitting. For simplicity, we have set equal to zero the 
imaginary parts of the NR cubic nonlinear coefficients in the 
appropriate resonant AEs. 

The FP solution of ( 20 ), i.e. with A m = 0, m - — , 0, +; 4> = 
0 with all 3 amplitudes constant and nonzero, when it exists, 
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Fig. 2. 2 : 1 Resonance: Amplitude modulations, as a function of the resonance detuning parameter <5^ in units of Ko; the large and small amplitudes 
correspond to modes a and b, resp. 


corresponds to a frequency-locked triplet. In addition to such a 
triple mode (TM) there can also be single mode (SM) and double 
mode (DM) FPs. Finally we recall that the resonant AEs can 
also have time-dependent solutions corresponding to amplitude 
modulated pulsations (BGS). 

In this section we explore what happens to a TM (frequency 
locked) FP as |<5 ^/ko| increases. This is of interest in a given 
star because of a variation of the growth rates among the excited 
multiplets, or, when we consider a collection of stars, because 
of the variety of rotation rates. 

As far as the corresponding NR AEs are concerned, there is 
no need to write them down, because they are simply obtained 
by setting the r m ’s equal to zero in Eqs. 20a, b and c. It is easy 
to show that the q and q are identical for this type of resonance 
(this can readily be shown following Buchler & Kovacs 1986a). 
Depending on the values of the parameters, the NR FPs in gen- 
eral can also come as SM, DM, TM, and may exist or not, be 
stable or unstable, depending on the coefficients. It is important 
to note the absence of the fourth equation for the phase, i.e. 
one no longer has an equation i = 0 as in the resonant case. 
This has as an important consequence that the NR triplet is 
necessarily unlocked. 

As in the previous section one can show that the frequency 
locked FP solutions cannot exist when is sufficiently 

large. What happens to these frequency locked solutions away 
from the resonance center? 

It should be obvious that a general parameter study, neither 
analytical nor numerical, is not possible and that we have to 
resort to numerical methods to provide illustrations of possible 
behavior. 

First example 

In a first such illustrative example we set the phases 6 m equal 
to zero for simplicity. Further, the coefficients in the AEs are 
chosen so that a stable resonant TM exists at 6u? = 0, on the one 
hand, and so that the NR AES have a stable TM as well. 

Our parameter values are: <7n = 0.326, q \ 2 = 0.2062, <713 = 
0.19, q 2 \ = 0.2122, <722 = 0.3425, q 2 3 = 0.2139, <731 = 0.2351, 
<732 = 0.2132, <733 = 0.3194, r_ = 0.51, r 0 = 0.018, r + = 0.6. 
For the growth rates: k_ = (1 + 6. x 10 -3 ) and k+ = (1 —5.x 
10“ 2 ). We have again rescaled the time in the AEs (Eqs. 20), 
by setting k 0 = 1 . 
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Fig. 3. Amplitudes and phase <£ of a resonant triplet ( = 1 as the 
resonance detuning parameter <?w/ko is increased; solid line: stable, 
dashed: unstable solution; Dotted lines: amplitudes of the NR triplet 
FP. 

In this simplified case, it can easily be checked that the res- 
onant TM FP satisfies the symmetry relation 4> — > — * 

— Su;, so that we can limit ourselves to positive 

With the chosen parameter values, one has a stable TM 
(frequency lock) right at the resonance center (Sjj = 0) with 
A- = 1.524 , Aq = 0.651, A + = 1 .367, 4> = 0.00. 

On the other hand, the NR AEs, that apply for large 
have the following solutions. A TM = 1.343, Aq = 1.152, 
A+ = 0.873. (Note the different asymmetries in the amplitudes 
of the resonant TM and NR TM FPs). The question arises how 
the system evolves from the resonant to the NR TM as <*w/ko 
is increased. 

Fig. 3 shows the results. The resonant TM FP remains stable 
(solid lines) as long as Suj/kq < 1 .367 above which it becomes 
unstable (dashed lines). For comparison, the amplitudes of the 
NRTM is indicated as horizontal, dotted lines. For large 6u>/kq 
the story is similar to the one described for the 2: 1 resonance 
in the previous section, namely the amplitudes oscillate around 
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the NR TM values with a period ~ 2i t/6uj. The magnitudes 
of the amplitude modulations decrease rapidly with [<5 u;/ko|- 
Ultimately, for large Suj/kq, the magnitudes of the modulation 
are so small that the star would be seen in practice to pulsate 
with constant amplitudes with the latter given by the NR AEs. 

What happens in the vicinity of the transition 8u)/ko 
1 .367 does not seem generic, and depends on the chosen pa- 
rameter values. In this example one obtains first a Hopf bifur- 
cation with oscillations about the now unstable resonant TM, 
but further away the oscillations occur about the NR TM, as 
just mentioned. The reason for this smooth transition from the 
constant amplitude, frequency locked pulsations to amplitude 
and frequency modulated pulsations seems to be that the reso- 
nant TM FP amplitudes nearly coincide with the NR TM ones at 
\6w/kq\ ~ 1.6 just after the onset of instability of the resonant 
TM solution. 

The transition from constant amplitude and frequency 
locked pulsations to modulated pulsations can be confirmed 


with a simple system of 3 coupled pulsators i.e. 


x+ 

= 2k—x + 8 R~yyz - 891 \±x 2 



-4<7i2(j/yz + y 2 i) - 4q l }(zzx + z 2 x) 

(21a) 


u>ly = 2 K 0 y + AR^y{xz + xz) - Aq 2 \(yx 2 + \ jxx) 



-2 q 2 jyy 2 - Mi^yz 1 + yzz) 

(21b) 

2+ 

uj 2 z = 2k+z + 8 R+yyx - 4q it (xxz + x 2 z ) 



-4qi 2 (yyz + y 2 z) - Sqxzz 2 

(21c) 


which have been chosen so that the amplitude equation for- 
malism applied to the system 21 yields the AEs of Eq. 20. 

To mimic for instance a rotationally induced spli tting, th e 
side frequencies are chosen as: u± = tuoO ± 0.5^6uj/u;o + 
0.5<5*j/u;o) where <5 u//lco oc (fi/u>) 2 . The same values as above 
are taken for the nonlinear coefficients q and r. The time in (2 1 ) 
has been scaled with u>o = T For our numerical experiments, 
two values of kq/u)q are considered below, viz. 10 -2 and lO -3 ,.. 

Fig. 4a shows the behavior of the component 2 with time 
for 6u)/k 0 = 0.5, 4.0 and 50.0. For better visualization we have 
taken kq/ljq = 1.5 x 10 -2 on the top. Transients have been 
removed and the time origin is arbitrary. At &j/ko = 0.5, 2 
pulsates with a frequency 2?+ with a constant (resonant) ampli- 
tude A+. At Sui/kq = 50, 2 pulsates with an almost constant NR 
amplitude. At 6uj/kq = 4, the pulsation is amplitude modulated 
with a modulation period of In/Su). 

The bottom figures show the Fourier amplitude spectra of 
the signal x + y + 2 built over a time interval of 240 basic periods 
P 0 - 2tt/lJq with 6uj/ko = 0.5. Here Kq/vo = 10 -3 . 

A Fourier spectrum of this signal (Fig. 4b) shows how the 
triplet goes from a frequency-locked to a frequency-unlocked 
pattern when the ratio of the departure from exact resonance (i.e. 
from equal splitting) to the growth rate increases, everything 
else kept constant. Note that at 6uj/k 0 = 4 in the intermediate 
regime, the time interval of the signal in Fig. 4 corresponds to 
about 1 period of modulation. It can be easily checked that for a 
longer time interval i.e. covering several periods of modulation, 
additional side frequencies appear in the Fourier spectrum. 



Fig. 4. top from left to right: Temporal behavior of the component 
z(t) (Eq.21c) for = 0.5,4,50 with ko/u>o = 10 - \ bottom: 

Corresponding Fourier amplitude spectra, see text. 

As far as observations are concerned, with typically 
Sut/ujQ ~ (ft/^o) 2 ~ I0 -6 , the rotational frequency, and a 
basic pulsation P 0 = 2tt/u>o — 200s, the period of modulation 
fora variable white dw’arf is of the order P mo d ~ 10 6 Po ^6-7 
years. For a delta Scuti star, with typically ~ 10 -4 and 

P 0 ~ 90 mn, P mo d ~ 10 4 P 0 ~ 1 - 2 years. Therefore, with 
the lengths of current continuous observation runs, if the the 
star were in the intermediate regime, only short phases of the 
modulation would be captured and the amplitude would barely 
vary over a same run. Of course, one would observe a pattern 
in amplitude which would vary from one run to depending at 
which phase of the modulation they are recorded. 

Second example 

Turning to the full case, that is with complex nonlinear co- 
efficients entering the AEs (BGS), leads to similar conclusions. 
The set of possible types of behavior for the full system, how- 
ever, can be more complex. The symmetry —> 

-6uj no longer exist and the system can behave quite differ- 
ently for and — In addition, several types of solution can 
coexist. 

Fig. 5 shows the evolution of the amplitudes of the resonant 
TM solution as one moves away from the exact resonance for 
the set of parameters 

q n = (0.326 — z 1.413); q ]2 = (0.2062 - i 1 .394); 
g 13 = (0.19 — i 3.368); q 2 1 = (0.2122 - i 1.351); 

<722 = (0.3425 - i 2.521); q 2 3 = (0.2139 - i 1 .412); 

<731 = (0.2351 - i 3.26); q 32 = (0.2132 - i 1.419); 
q 33 = (0.3194 -i 1.451); 
r_ =(-0.018 - *0.51) = 0.5 103exp( 1.5355 0; 
r 0 = (0.024 - i 1 .05) = 1 .0503exp(- 1 .5479 i); 

7> =(-0.021 - *0.6) = 0.6004exp( 1.5358*)- 
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Fig. 5. Amplitudes and phase of a resonant triplet I = 1 as the ratio 
1 / /vy | increases - example 2; (solid line: stable and dashed lines : 
unstable triplet); dotted line : amplitudes of the constant amplitude, 
nonresonant triplet which exists and is stable far away from the reso- 
nance. 


In the resonance center (Su: = 0) the resonant AEs have a 
TM (frequency locked) FP. Here the solutions and the stability 
diagram are no longer symmetric about the resonance center. As 
in the previous two examples the frequency locked FP becomes 
unstable for a sufficiently large deviation from resonance, i.e. 
for 6u>/kq < -9.033 and 6u>/kq > 2.53, respectively, as Fig. 5 
shows. For these numerical coefficients the system bifurcates 
abruptly from the frequency locked FP to an amplitude mod- 
ulated pulsation. The amplitude modulations occur about the 
NRT FP solution, at first with large excursions which later scale 
as \/6uj and whose period scales also as 1 /6u>. However, as 
for the 2:1 resonance example above, there is hysteresis, i.e. a 
narrow regime of 6u>/kq in which the frequency locked and the 
amplitude modulated solutions coexist. 


4. Resonant multiples 

We assume now that two excited multiplets can interact res- 
onantly, but that they are decoupled from other possible ex- 
cited multiplets. In the following we label the two interacting 
multiplets a = {a m , m = — 4?”' >4} an d b = {6 m , m = 

- 4 , 

Among all possible two-mode resonances, three of the types 
we consider below are especially important because they can 
affect the linear stability of a pulsation. These are: integer 
resonances (Sect. 3.1), direct resonances (Sect. 3.2), and half- 
integer resonances (Sect. 3.3). 


4. 1. Integer resonances: ncj a « u>b 

For didactic purposes we first consider the simplest integer res- 
onance; that is, the 2: 1 resonance. For simplicity of notation we 
will temporarily use the symbols a and b to stand generically 
for any component of the associated multiple!. 

** 2:1 resonance - 2uj a ~ Ub 

The appropriate AEs are easily seen to be (e.g., Dziem- 
bowski & Kovdcs 1984, BK86a) 

a = o a a - q aa aa* a - q a bbb*a, + r\ba * 
b = &bb - qba aa *b - qt,bbb*b + r 2 a 2 

where an asterisk denotes complex conjugate. The four cubic 
q a $ coupling terms are the important amplitude saturation terms 
that always arise, whether there is a resonance or not and whose 
expressions are given in BG84 and GB94. In all our applications 
of AEs to radial pulsators we have found that the real parts of the 
q a p coefficients are always positive (The physical reasons are 
addressed in Buchler 1 993). In some of the following discussion 
we shall assume, for specificity, that this also holds true for 
nonradial pulsators. 

The expressions for the resonant coupling coefficients r\ 
and r 2 are given by the integrals over space and solid angle 

r\ oc rj oc (a]N 2 |a*b) m 

oc / drpr 2 f (IQ C(r)Yl N 2 £»&(r) YJ. Y ( , . J 

where Y£ is the usual spherical harmonic and Nj is an operator 
acting on £ (the r-dependent part of the eigenfunctions) which 
comes from the second order Taylor expansion of the hydrody- 
namical stellar equations (BG84; GB94, and see Sect. 5). Be- 
cause of the assumed slow rotation and accompanying spherical 
symmetry of the equilibrium stellar model, the stellar structure 
operators 2T (the linear pulsation operator), N 2 , N3, • • ■, have 
even parity and azimuthal symmetry. The angular integrals thus 
restrict the types of couplings that are allowed (Dziembowski 
1984, GB94). Therefore the angular part of (23) simplifies to 

r, oc r 2 * « [ YlYtY th dSl. (24) 

Parity considerations (see Sect. 5) show that the coupling terms 
r\ and ri vanish identically 

(a) when 4 is odd and W a , or 

(b) when 4=0 and 4 ¥ 0. 

Of course, additional constraints on the m a and nib values arise 
from angular momentum ( L z ) constraints. Here, = 2m a 
must be satisfied. 

It is of interest to note the consequences of the respective 
linear stability or instability of the two modes. If mode a is lin- 
early unstable ( K a > 0) and a 4 0, then mode b will necessarily 
be nonzero, even when it is linearly stable. (This is the situation 
that occurs in the bump Cepheids, BK86.) In other words, mode 
a always entrains mode b . 

When mode b is unstable with B ¥ 0, then A = 0 is always 
a solution. One can easily show that when this FP solution is 
unstable then a solution with .4/0 and B ^ 0 exists. In that 
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Table la. Resonance 2u> Q ~ 


L 

h 

allowed 

allowed (m fl ,m&) 

0 

1 

no 


0 

2 

no 


1 

0 

yes 

(0,0) 

2 

0 

yes 

(0,0) 

I 

1 

no 


2 

1 

no 


1 

2 

yes 

(0,0); (±1, ±2) 

2 

2 

yes 

(0,0); (±1,±2) 


Table lb. Resonance 3u> a ~ 


la 

ib 

allowed 

allowed (m a ,mt>) 

0 

1 

no 


0 

2 

no 


1 

0 

no 


2 

0 

yes 

(0,0) 

1 

1 

yes 

(0,0) 

2 

1 

no 


1 

2 

no 


2 

2 

yes 

(0,0) 


case one says that mode b parametrically excites mode a (cf. 
below under 1 :2 resonance). 

In Table la we have summarized the situation for the com- 
binations of the lowest 4 and 4 values that are of practical 
interest. 

In the realm of low £ nonradial modes condition (a) shows 
that the higher frequency mode thus needs to be an £ 5 = 2 
(even) mode, but that the lower frequency one could be £ a = 1 
or 4 = 2 . 

Condition (b) shows that a radial mode cannot excite 
a nonradial mode through this type of resonance unless 
the radial mode is the high frequency partner in the resonance. 
Thus, for variable white dw’arfs for instance, g-modes might be 
excited this way, but never p-modes. 

** n:l resonance - nuj a « n > 1 
We now go on to the general case 


a = 

b = 


+ r\ba* n ~ ] 
+ r 2 a n . 


£(a|N 2 |/3fc)-F(fc|N 2 | 7 5) 


where the denominators D k represent sums and differences of 
linear frequencies (e.g. BG84, GB94, BGS). Typically these in- 
direct terms are of the same size as the direct term. For n - 3 
one can still obtain relatively simple explicit expressions for 
the cubic resonant coefficients (as for the 1:1 direct resonance, 
BGS). For higher order (n > 3) resonances, however, the indi- 
rect terms are quite complicated. 

But, fortunately, we can take advantage of a general result, 
established in Sect. 5. 

When the direct coupling term vanishes because of par- 
ity and angular momentum (L z ) considerations then each 
of the possible indirect coupling terms also vanishes , and 
conversely . 

As a consequence, whatever the resonance is it is sufficient 
to study the simple direct term to determine whether a 
specific resonant term should appear in the AEs or not . 

This enables us to obtain the allowed n: 1 resonant couplings 
as listed in Table lb for the example of a 3: 1 resonance and for 
couplings involving / = 0 to I = 2 modes. 

The above constraints on £ values given in Tables la and 
lb show, for example, that the Blazhko effect in RR Lyrae stars 
cannot be explained as the coupling of a radial mode with 
a nonradial p-mode of higher frequency through a 2:1 or 3:1 
resonance nor, for that matter, through a rc:l (n > 1 ) resonance 
in general. 


4.2. Direct or 1:1 resonance: 




! ^6 


A special situation arises for a direct 1:1 resonance which, like 
the 3: 1 resonance, also contributes cubic terms. 

Denoting the amplitudes again generically by a and b we 
obtain the AEs 

a = . . , + r\b 2 a* + s\baa* + t\a 2 b* + U\ bbb* 

b = . . . + r 2 a 2 b* + s 2 abb * + t 2 b 2 a * + u 2 aaa* 


(27) 


(25) 


Here there are four possible resonant terms in each AE. As we 
have just seen, on the basis of the theorem of Sect. 5, we only 
need to consider the direct terms to determine which of these res- 
onant terms survive parity and angular momentum constraints. 
For these direct terms not to vanish, the constraints on I and m 
are the same as those given by the product of the four appropri- 
ate spherical harmonics. Thus the conditions for the existence 
of r\ and r| can be obtained from 


where, here and from now on, we omit the usual linear and cubic 
saturation terms. 

The angular integrals show that the “direct” coupling terms 
(due to N n ) vanish when 

(a) n£ a + 4 is odd, 

(b) 4=0 and 4 t^O. 

For n > 2, there are also “indirect” contributions to r x and 
r 2 . For example, in the 3: 1 case where the direct contribution is 
(a|N 3 |/? 7 < 5 ), the indirect contributions all contain terms of the 
general form 


/ 


Y £ * Y! ■ Y £b Y ib dn. 


(28) 


This implies that r x and are different from zero independently 
of the values of 4 and 4 - 
Similarly, 


/ 


Y? a Y /a Y 4 .Y IdQ 


(29) 


implies that sj, tj , and it 2 vanish when l a + 1 b is odd, or when 
( a = 0, V4. Finally, the integral 


(26) / VYe b Ye h Yl<m 


(30) 


170 


J.R. Buchler et al.: On the role of resonances in nonradial pulsators 


Table 2. Direct Resonance u ; a ~ 


(a 

b 

n* rz 

5 1 , t\, it 2 

U 1, $2, ti 

0 

1 

yes 

no 

no 

0 

2 

yes 

no 

yes 

1 

0 

yes 

no 

no 

1 

1 

yes 

yes 

yes 

1 

2 

yes 

no 

no 

2 

0 

yes 

yes 

no 

2 

1 

yes 

no 

no 

2 

2 

yes 

yes 

yes 


Table 3, Resonance (2n+l)u/ a « 2 a>6, n > 0 (e.g. 3u; a « 2wf>, 

5u/’a ~ 2uJ b ) 


e. 

tb 

allowed 

0 

1 

yes 

0 

2 

yes 

1 

0 

no 

I 

1 

no 

1 

2 

no 

2 

0 

yes 

2 

I 

yes 

2 

2 

yes 


The 2:(2nfl) case is obtained by inverting the roles of a and b 


yields the conditions that u\ , si> and f 2 vanish when either 4 +4 
is odd, or when 4 = 0, V4- 

The m-constraint is the same for all coefficients; i.e., m a = 
mb 

In Table 2 we examine the combinations of low £ 1:1 mul- 
tiplets of potential interest in nonradial stellar pulsators. 

If two multiplets are of the same type and one is a £ - 1 
and the other £ = 2, then Table 2 shows that all the 5, t and u 
vanish and there remains only one type of resonant term. This 
same simplification occurs when an £ = 1 p- or g-mode is in 
direct resonance with a radial mode. 


4.3 . Half-integer resonances: (2n+l)u; a ^ 2 u;*, 

One notes that this simplest 1:2 {n - 0) case is essentially the 
same as a 2:1 resonance provided that we invert the roles of a 
and b . The coupling terms thus vanish when (a) 4 is odd and 
V4, or (b) when 4 = 0 and 4^0. 

The higher frequency mode (a) must therefore have even 
4 to be able to parametrically excite the linearly stable lower 
frequency mode 6 (which can have any 4 value). Generally, a 
radial mode could thus excite any resonant g-mode, even when 
the latter is linearly stable. This situation is very similar to the 
three-mode resonance uj a % u ;& + u) c considered by Dziem- 
bowski (1982). 

The relevant terms in the AEs are 


d = * - + rifr 2 a* 2n 

b = • « * + r 2 a 2n+1 6*. 


(31) 


The direct contribution to r\ and r 2 is easily seen to vanish when 
4 is odd. This resonance can couple only central modes of the 
multiplets (m a = = 0). 

Eqs. 24 show that depending on the linear stability and the 
values of the nonlinear coupling coefficients three types of so- 
lutions may exist. 

When mode a is linearly unstable ( k q > 0), but mode b is 
linearly stable, and if furthermore the FP solution (.4 0. B = 

0) would be stable in the absence of the resonance (i.e., if + 
Qba^ 2 < 0), *hen mode b gets excited parametrically when (cf. 
Moskalik & Buchler 1990) the Floquet exponent 

Kb - q*A 2 + J\r 2 \iAM"+ 1) - ($u/ - C 42 ) 2 (32) 

is positive. The superscripts R and 7 refer to real and imaginary 
parts, respectively, of the saturation coefficients. If, however, 
mode b is also linearly unstable, then of course a double-mode 
solution occurs. 

On the other hand, when mode b is linearly unstable, but 
mode a is linearly stable, one easily sees from Eq. 3 1 that mode 
b cannot excite mode a. 

In Table 3 we summarize the situation for the (2?v + 1):2, 
n > 0, resonance. One notes that a radial a mode can couple to 
and parametrically excite a higher frequency 4 = 1 or 4 = 2 
mode, in contrast to the n:\ resonance. 


4.4. Other potentially important resonance, 4:3 4 jj a ^ 3*^ 
The AEs in this case contain the resonant terms of the form 


a = ... + rib 1 ' a * 3 

b = ... + r2<2 4 &* 2 . 


(33) 


Angular momentum constraints make the direct contribution 
to the r coefficients vanish only when mode b has an odd £. 
In addition, the coupling is possible only between modes with 
m = 0. Parametric entrainment of one of the modes by the other 
is not possible here because this resonance does not affect the 
linear stability of either single mode FP. It is only when both 
modes are excited that the resonance term is allowed to play a 
role. 


4.5. Three -multiplex couplings 

AEs can also be useful in the case of three-mode or multiplet 
couplings. Dziembowski (1982), for example, considered the 
coupling of three modes through a resonance of the type ^ 
u?b +u; c , which is well known in plasma physics under the name 
of resonant wave interaction. In later work Dziembowski & 
Krolikowska (1985) suggested that it is this type of resonance 
that limits the amplitudes of 6 Scuti type stars to their relatively 
small observed values. 

Also we have already noted that the situation of a rotationally 
split multiplet is mathematically equivalent to a resonance of the 
type 2u) a « uJb + u; c . 
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5. Properties of the coupling coefficients 

This section is very technical and can be skipped in a first read- 
ing. We prove here a very useful property of the nonlinear 
coupling coefficients. We assume that the equilibrium model 
is spherically symmetric (slow rotation, no magnetic fields). In 
zeroth order in the rotation rate O the linear eigenvectors are 
then represented by spherical harmonics Yi m which have parity 
(-)*; i.e., they obey II Y im = and L z Y £m = mY em 

(e.g., Blatt & Weisskopf 1952). Because of the spherical sym- 
metry the nonlinear structure operators have positive parity as 
well as m = 0. 

To briefly introduce what follows, we first write down the 
pulsation equations in a compact form using the convenient 
Dirac notion where z represents, for example, the density or 
other stellar variable in an Eulerian representation (as in BG84 
or GB94): 

^ = .Z|z)+N 2 [|z)] + N 3 [|z)] + --- (34) 

where Z is the linear pulsation operator and the operators N 2 , 
N 3 , etc., are quadratic, cubic, etc., in the displacement vector z. 

The AE formalism then provides a sytematic expansion of 
z as 

|z(t)) = |z 0 (t» + |zi(t)} + |z 2 (t)) + ■•• (35) 

where 

Mt)> = 5^aae ,iJo< |a) (36) 

a 

and |zj(t)), |z 2 (t», etc., are quadratic, cubic, etc., in the am- 
plitudes {a a (f)}, and thus represent modal interaction terms 
(“combination frequencies”). The sum is over all dominant 
modes q. 

The resonant coupling coefficients, which measure the 
strength of the coupling between modes, are given by such inte- 
grals as in Eq. 23 and Eq. 26 and are broken up into “direct” and 
“indirect” contributions, as indicated by the discussion leading 
to Eq. 26. 

Theorem : If, genetically, the direct coupling term is zero 
because of parity and angular moment um considerations , 
this implies that each of the possible indirect terms is also 
zero , and conversely . 

Parity. Consider first the lowest — that is, the cubic coupling 
term — for which there is a direct and indirect contribution. The 
direct term is of the form 

r dir oc (a|N 3 |bcd). (N 3 zJ) (37) 

The expression in parentheses indicates the perturbative origin 
of this term; cf. expansions (34) and (35). The parity of the 
RHS is Il(bcd) = n(b)n(c)IT(d), and the parity of the stellar 
structure operator is II(N3)=1 . The direct term thus vanishes un- 
less 11(a) = n(c)II(c)n(d). The indirect term (Eq. 26) contains 
terms of the form (e.g. BG84, GB94) 

r ind * Y <a|N 2 |bk)<k|N 3 |cd>. (N 2 zo Zl ) (38) 

k 


For the purpose of clarity we have ignored denominators (sums 
and differences of frequencies) and factors of unity. The sum 
runs over all eigenstates of the system. This is one of the rea- 
sons why it is so difficult to make an ab initio computation 
of these terms. It is understood that b, c and d represent both 
bed and their complex conjugates, and that the appropriate per- 
mutations of (bed) are included. The last factor in (38) van- 
ishes unless n(k) = TI(c)n(d), and the first one vanishes unless 
11(a) = n(b)II(k) = n(b)II(c)n(d); i.e., the same condition as 
for the direct term. 

Consider now the quartic coefficient. The direct term 
rdir cx (a|N 4 |bcde) (^zjJ) (39) 

vanishes unless 11 (a) = n(b)TT(c)n(d)n(e). 

The indirect terms are more complicated now. They are 
given by 

r X nd oc y (a|N 3 |bck)(k|N 3 |de) (N,z*z,) 

k 

+ (a|N 2 |bk)(k|N 3 ]cde) (N 2 z<>z 2 ) 

+ YJ2 ( a l N * l bk > < k i N 2 l ck '> ( k ' INj |de) . 

k k' 

(N 2 z 0 z 2 ) (40) 

It is easy to see again that all three contributions vanish unless 
11(a) = rT(b)n(c)IT(d)n(e). In fact one sees that each nonvan- 
ishing matrix element successively “passes through” the parity 
from the left side to the right side. This holds true for all matrix 
elements, whatever their order. Since all higher order indirect 
terms are built up in a manner similar to Eq. 40 we conclude that 
the parity part of the theorem is valid for all orders of coupling. 
Angular momentum: The nonvanishing of a given matrix ele- 
ment requires that the L z component on the right match the L z 
component on the left; i.e., mjhs = m rhs (because the all 
have 7Ti - 0). 

Referring to the previous section on parity, mutatis mu- 
tandis we see that this sum is passed through the nonvanishing 
matrix elements. Each of the indirect contributions therefore 
vanishes generically when the direct term vanishes and vice 
versa. QED. 

The theorem thus implies that if one is interested in deter- 
mining which of the normal coupling terms survive in the AEs 
it is sufficient to consider the expression of the simple, direct 
term. 

6 . Possible resonances in observed power spectra 

In this section we briefly review some examples taken from the 
observational literature which, in many cases, strongly suggest 
that resonant interactions are indeed present in nonradially pul- 
sating stars. For the most compelling examples we also outline 
some consequences of the suspected interaction between modes. 

For some of our examples we must caution the reader that 
the quoted observational results depend on data that may be mis- 
leading because of insufficient time coverage of the star being 
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discussed. For these stars, and others of course, this section then 
serves to draw attention on the one hand to the need for and on 
the other to the potential benefits of more extensive observations 
directed to the study of these resonances. 

6.1. Integer resonances: ruv Q « iCt 

The ZZ Ceti star GD 385 (Kepler 1984) provides a perfect ex- 
ample that can be explained by a 2:1 resonance. The Fourier 
spectrum of GD 385 has three well-defined peaks located at 
/i _=3,9043, /i+=3.901 2 and / 2 =7.8055 mHz. The observations 
show that 

h =/!-+/.♦ < 41 > 

to a relative accuracy of 10 -6 . In the following we use the sym- 
bol / for observed frequencies which are to be interpreted as 
nonlinear frequencies (/ = u>/ (27 t)). 

Let us assume, as suggested by Kepler (1984), that the two 
low-frequency peaks are modes ±m belonging to the same mul- 
tiple! of given 7 and radial order. The central component of this 
multiplet is frequency locked to / 2 to within one part in 10 -6 . 
The linear eigenfrequencies of the two interacting multiplets 
then are in a near 2: 1 resonance with u) 2 ~ 2cc\o ~ ctZ] _ + uj\+. 
The 777 = 0 central peak f\o is not seen, perhaps due to inclina- 
tion effects, or due to nonlinear effects: the m = 0 mode indeed 
may or may not be involved in the pulsation. Table la shows 
that the higher frequency multiplet containing / 2 must therefore 
be even and probably i 2 “ 2. (An l of three or higher has never 
yet been identified in a variable white dwarf) and m 2 = 0. The 
lower frequency doublet, on the other hand, can be of odd or 
even angular momentum and the likely possibilities are £\ = 1 
or 2. In either case, |roj | must be equal to 1 . 

One might think that / 2 is just a simple harmonic of / ]0 , with 
the /to peak presumably not observed because the inclination 
angle is close to tt/ 2. We believe this to be unlikely because the 
harmonics of f\± should then have similar pow f er, but they are 
not visible. 

One can easily show that it is the phase $ = (f>2 — <P\- ~ 
0i+ which enters the apposite AEs. If the pulsations correspond 
to fixed points of the AEs then, as we have seen in Sect. 2, 
this leads to the nonlinear frequency lock of Eq. 41. However, 
some amplitude variations have been reported by Kepler (1984) 
that occur on the very long time scale of years. Would this be 
confirmed, it would imply that it is not a fixed point solution, 
but rather a time-dependent one. Such time-dependent behavior, 
if generated by the resonance, would of course not vitiate the 
angular momentum attributions of the mode. 

But because the observed modulation time scale is very long, 
the pulsations may still appear to be in a frequency lock when 
observed over shorter scales. Indeed, as the AEs show (Sect. 3), 
depending on the values of the coefficients, a resonant coupling 
can generate modulations of amplitude that are quite different 
in size from the period (or phase) modulations. In the present 
case, the magnitude of the amplitude modulation would be much 
larger than that of the phase modulation. 


6.2. Direct resonances: uj q ^ 

During the course of evolution of a star, pulsation frequencies 
change and, in some circumstances, the frequency of one mode 
may closely approach that of another mode of the same t. In this 
case, called avoided crossing, the two linear eigenfrequencies 
never become exactly equal to each other. Near the resonance, 
the two modes take on a mixed character having some of the 
properties of both modes. This avoided crossing, which is a 
purely linear effect, is discussed in a stellar context by, for ex- 
ample, Unno et al. (1989, Sect. 15.3.1). Nonlinear effects can 
however cause the two frequencies to lock, i.e. for the tw'o fre- 
quencies to coincide exactly. 

Possible examples of direct resonances may be present in 
the variable white dwarf PG 1 1 59-035 where only high-order g- 
modes of P = 1, 2 are observed. With a rich spectrum of modes 
(125 were resolved in PG1159 by Winget et al. 1991) there 
may be an occasional overlap in frequency between an ( = 1 
and an ( = 2 in the Fourier spectrum. This is especially true at 
low frequencies where, for g-modes, the frequency spectrum 
becomes particularly dense. Examples of overlap are reported 
in Winget et al. (1991, Table 3) where an P value is assigned as 
either P = 1 or 2. Reference to our Table 2 shows that a direct 
resonance is indeed possible here (with only the r\ and r 2 terms 
contributing) which may distort what is normally expected to 
be seen in the two overlapping multiplets. To disentangle the 
effects here is both an observational and theoretical challenge. 

Again, depending on the physical situation, and thus on the 
values of the coupling coefficients, a direct resonance can lead 
to unsteady behavior (amplitude and phase modulations). If 
the modulation amplitude is small and sinusoidal, a single pul- 
sation then appears in a power spectrum as if it were governed 
by three or more closely spaced frequencies. 

Finally we note that a direct resonance between a radial and a 
nonradial mode has recently been suggested as a possible cause 
of the Blazhko effect in RR Lyrae stars (Kovacs 1994, 1995; 
van Hoolst 1995). 

6.3. Half-integer resonances: (27? + I)u> a « 2c ot 
6.3.1. Variable white dwarfs 

PG 13513-489 

The light curve of the DBV star PG 135 1+489 described by 
Winget et al. (1987) is dominated by a single period of 489.5 s 
corresponding to a frequency of /o « 2.05 mHz. Strong har- 
monics are also seen at 2/o, 3/o, etc. Of interest to us, however, 
are lower amplitude signals seen in most of the observations 
made in the late 1980’s at (approximately) 1.47/ 0 , 2.47/ 0 , and 
3.47/o, in order of increasing frequency and decreasing ampli- 
tude. Goupil et al. (1988) have claimed that power is also present 
in the data at 0.47/ o in two observations reported by Winget et 
al. (1987). The situation is murky, however, because red noise 
due to low-frequency perturbations in the atmosphere can mimic 
or even cancel out the signal we seek. This is especially serious 
if, as in the Winget et al. (1987) observations, data are taken 
from a single site with the inevitable one day, or longer, gaps. 
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To circumvent this, as yet unreported observations utilizing the 
“Whole Earth Telescope” (WET) were made in early 1995. (For 
a description of WET seeNatheret al. 1990.) With the extended 
global coverage capabilities and improved resolution of WET, 
a signal at 0.46/o is seen very clearly as are ones at 1.46/o, 
etc. The amplitude of the peak at 0.46 /o is only 4% of / 0 itself, 
which illustrates the virtue of extended coverage in resolving 
peaks of low amplitude. Although the 1995 data have not yet 
been completely reduced and analyzed, the peak at /o shows 
evidence for splitting (or unsteadiness) whereas, to the limits of 
the observations, neither 1 .46/o nor 2.46/o seem to have multi- 
plet structure. (There are other low amplitude signals in the data 
but they appear to be independent modes not directly associated 
in a resonant fashion with the main peaks.) 

It seems reasonable to interpret the Fourier spectrum as 
generated two modes, at /o and at f\ = 1.46/o (or perhaps 
/i = 2.46/ 0 ). All other peaks would then be harmonics of / 0 
and combination peaks of /o and /] . 

There is no absolutely compelling reason why these could 
not be two independent, nonresonant modes. It is difficult to 
understand though why there would be just two excited modes, 
with a rather tight frequency ratio of 1 .46. A coincidence seems 
unlikely because there are two additional white dw ; arfs known 
with a similar Fourier spectrum and a frequency ratios of 1.52 
and 1 .54, respectively (see below). However, the observed be- 
havior appears quite natural if we assume that a half integer 
resonance is active in this star (3:2, or perhaps 5:2). 

It might be objected that the observed intermediate frequen- 
cies do not appear at 1 .5 (or 2.5), but that they always fall short 
of those values. In a resonance situation one would expect non- 
linear effects to lead to either of tw f o situations, viz. either exact 
frequency-lock, or fluctuations about the resonance condition 
i.e. in both directions. The spectrum of PG 1351 can however be 
explained if the observed second frequency (f\ =1.46 or 2.46 
/o) is one of the m ^ 0 components of the multiples and is the 
only one that is large enough to be observed. The work of BGS 
has shown that large amplitude asymmetries within a given mul- 
tiple! can easily happen. In the present case, this would yield a 
rotational splitting ± ~ 60 — 80 fiHz of the I.5/o peak. This 
indicates a rotation that is probably too fast if the f\ peak cor- 
responds to the 777 = 1 component (part of the split could be 
explained by nonlinear effects). The size of the split therefore 
favors again an £ = 2,m = 2 (or higher £ y m) component. In 
either case, if we are dealing strictly with nonradial modes of 
low order — and the observed periods do rule out £ - 0 modes 
- then Table 3 suggests that the mode at / 0 is an £ - 2 mode. 
(Note that £ > 3 modes have not yet been identified in white 
dwarf light curves.) 

With either of the two resonant possibilities mentioned 
above, Table 3 restricts the mode at /o to be an £ = 2 mode.. 
(The observed periods definitely rule out £ = 0). The discussion 
surrounding Table 3 also states that a 3:2 or 5:2 resonance can 
couple only m - 0 modes of a possible multiplet. 

If /o, with the aid of future observations, turns out to have 
£ = 2 then that would be a strong argument for our resonant 
identification. 


GD154 

Another curious example is the DA variable white dw arf 
GD 1 54 first reported on by Robinson et al. ( 1 978). This 1 978 pa- 
per describes a star that was a virtual tw in of PG 1351 discussed 
above except that the ratio of near 3:2 (or 5:2) was 1 .52/ 0 , etc., 
rather than (in the latest observations) 1.46/ 0 , etc. With this 
information alone, we would assign £ = 2 to the mode at /o 
(where, in this case, /o = 0.842 mHz). 

If, however, we only had the more recent WET observations 
of Pfeiffer et al. (1993) we w ? ould tell a different story: they 
observe no half-integer signals at all and we would have no 
reason to assign £ = 2 to / 0 . The true story, wc believe, is that 
the assignment £ = 2 may still be valid but — and being true to 
the whole story we tell here — nonlinear effects of the 3:2 (or 
5:2) resonance force the half-integer mode to essentially appear 
and disappear on time-scales that are as yet unknown. Pfeiffer 
et al. may just have caught the star at an inopportune time. 

BPM31594 

This is a white dwarf to have been observed with simi- 
lar Fourier spectral features (CFDonoghue, Warner & Cropper, 
1992). Here the second frequency appears at 1.54 / 0 , and har- 
monics at 2.54, 3.54. 

6.3.2. 6 Scuti stars 

Another potential candidate for a 3:2 resonance is afforded by 
the 6 Set star DL Eri (HR 1225) reported on by Poretti (1989). 
The Fourier spectrum seems to be dominated by three fre- 
quencies, /i=6.41 , /2_=8.98 and /2+=10.26 c/d. One notes that 
3/i = fi— + fi+ to within the claimed accuracy of the frequen- 
cies; in other words, the low frequency is in exact resonance 
with the center of the high frequency doublet. Table 3 indicates 
that the low frequency mode must have even £ (and could also be 
radial); there are no £ constraints on the high frequency mode. 
It remains to see whether these conclusions would survive an 
increase in resolution and signal-to-noise ratio as provided by a 
multisite observation campaign. 

6.4. 1:1:1 resonance: rotational splitting 
6.4. 1 . The ZZ Ceti star G226-29 

This star is observed to have a Fourier spectrum composed of 
a triplet of modes which is split almost exactly evenly in fre- 
quency (Kepler, Robinson & Nather 1983). No amplitude vari- 
ations are mentioned. More recent observations using WET by 
Kepler etal. (1995) provide updated information. The main peak 
appears at P = 2n/uj = 109.3s and has an average splitting of 
16.15^Hz. It appears that the frequencies in the triplet have re- 
mained constant to within about one part in 10 4 (as the limit 
of their measurements) and there is only marginal evidence that 
one of the amplitudes has changed by perhaps 20% . From these 
considerations it seems that the theoretical connection of am- 
plitude constancy and frequency locking is indeed verified for 
this star. 

Since white dw'arfs are slow rotators and since the linear 
deviation from equidistance is second order in the rotation fre- 
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quency (Eq. 1 7 and following paragraph) nonlinear effects need 
not be strong to bring about frequency locking. Thus, in this case, 
the splitting gives a good estimate of the rotation frequency (Ke- 
pler etal. 1995), from which we infer 6u>/uj % 1.2 x 10 _5 .The 
approximate condition for frequency locking (Eq. 17) requires 
a growth rate to frequency ratio k/u) of the order or larger than 
~ 10 -5 — 10 -6 . This figure is consistent with what is currently 
obtained for low overtone g - modes of DAV models (Bradley 
1994). 

6.4.2. DOV White dwarfs 

Two additional variable white dwarfs are of interest for this dis- 
cussion and these point to the need for more detailed studies 
of 1 : 1 : 1 resonances. These stars are PG 1 1 59-035 (Winget et al. 
1991) and PG2 1 3 1 +066 (Kawaler et al. 1 995), which have been 
studied extensively with WET. They are similar with respect 
to luminosity and effective temperature (luminous and hot) and 
both show conclusive evidence for a large number of t = 1 
nonradial triplet modes. In PG1 159 the triplet of highest am- 
plitude has undergone a modulation in amplitude (of unknown 
time scale) among members of the triplet, while some triplets 
in PG2131 have shown curious changes in relative frequencies 
within a triplet. In PG2131 the low frequency peaks appear 
sharp, whereas the central and especially the high frequency 
peaks show a broader structure. The observed broadening of 
the lines in the Fourier spectrum could be caused by amplitude 
modulations which can have different depths for the sub-peaks 
of a given multiplet (BGS). 

The Fourier spectra of PG 1 1 59 and PG2 131 are dominated 
by just a few modes, but without any apparent resonances. All 
the other peaks have much lower power. According to Kawaler 
et al. (1995) these low-amplitude modes are observed to come 
and go on time scales of months to years. This is the kind of 
time scale for amplitude modulations that one may expect on 
the basis of linear growth rates in many circumstances (Lee & 
Bradley 1993). Could these low-amplitude intermittent modes 
be stochastically excited by, for example, turbulence associated 
with convection? 

We mention the following as a possible scenario for the 
observed spectra. One or a few modes (multiplets) are linearly 
unstable and grow to saturation amplitudes given by 

A 2 » ?-(« ± x/k 2 -4q. c S) * - (42) 

2 q q 

where k and q are a typical growth rate and cubic saturation 
coefficient (see, e.g, Buchler, Goupil & Kovdcs 1993). 5^ is 
the noise intensity that drives the mode and has its origin in 
the convective and turbulent motions. The intra-multiplet am- 
plitudes primarily owe their variations to nonlinear effects that 
are inherent in the 1:1:1 resonance (BGS) of the multiplet with 
itself. 

The modes on either side are increasingly stable (k < 0). If 
they are stable but stochastically excited, their amplitudes are 
given by A 2 j=s f\n\. (This is the usual linear result valid 
when q< 5^ /k 2 <$: 1). Thus assuming that the noise intensity is 


more or less the same for neighboring multiplets, one expects 
the amplitudes to taper off gradually on either side of the linearly 
unstable modes. It was shown theoretically (cf. Fig. 1 of BGK) 
that the amplitude fluctuations are largest for modes near k - 0. 
If the above picture is correct then one expects the modes of 
lower amplitudes to undergo larger amplitude fluctuations than 
the dominant modes of large amplitude, with a concomitant 
correlation with broad and narrow peaks in a Fourier spectrum, 
features w'hich are amenable to observational verification. 

Superposed on this scenario would be the effects of trapping 
on the modes. Indeed one expects the growth rates to “oscillate” 
as a function of modal overtone. Thus there could be one or more 
further enhancements of the amplitudes on one or both sides of 
the dominant peaks. The spectrum of the DB variable white 
dwarf GD358 (Winget et al. 1994) shows amplitude variations 
between successive modes that might be an expression of this 
scenario. 

6.4.3. 8 Scuti stars 

The observed frequencies f\ = 7.217139, /> = 7.346146, fo - 
7.455268 c/d of the 8 Scuti star 1 Mon are locked, i.e f\ + fo = 
2 / 2 , to a relative accuracy of 10 -5 . The frequencies f\ and 
have been identified as the m = — 1 and m = +1 components 
of a £ = 1 triplet split by rotation. The rotation rate is then 
ft = 0.131 revolution per day (Balona and Stobie, 1980). That 
f 2 is a radial mode as it has been identified or the centroid mode 
of the / = 1 triplet is still questioned (Balona and Stobie, 1980; 
Breger, 1992). 

Whether the observed frequency-lock is in agreement with 
Eq. 1 7 depends on the eigenfrequencies and on the growth rates 
of the modes. If fi is the frequency of the m = 0 component of 
the triplet, the frequency lock ought to exist for a large enough 
growth rate i.e. k/jj > (ft/a ;) 2 ^ 3.210“ 4 ( u ; « 0.5(/j + 
/ 3 )), again within a factor 10 due to the involved, but discarded 
linear and nonlinear coefficients. This is the order of magnitude 
expected for the strongest linearly excited modes of 6 Scuti 
stars. This would explain why the triplet appears equally spaced 
when a linear framework predicts, with the above rotation rate, 
a significant departure from equidistance. 

On the other hand, the selection rules show that a locked 
radial mode, associated with fi , is also possible. Of course, this 
requires that the £ = 0 and the (unseen) m = 0, t = 1 modes 
are close enough for the resonance to be effective in locking the 
frequencies, a closeness that is encountered relatively often in 
models of 8 Scuti stars. On a practical side, this second possi- 
bility shows that an observed multiplet of locked frequencies 
is not necessarily associated with modes of same degree ( but 
may well correspond to modes of different degrees. 

7. Conclusion 

We have discussed the constraints imposed by resonances and 
have considered the observational material to determine where 
such resonances may play a role, and where new information 
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may be gained from applying the amplitude equation formalism 
to nonradially pulsating stars. 

We have assumed that individual or perhaps two or three 
multiplets can be studied in isolation from other observed mul- 
tiplets. We would like to point out that this assumption may be 
even better justified when the additional observed multiplets, 
rather than being excited because of positive driving (by the 
K-mechanism, for example), are merely stochastically excited. 

With these assumptions we have shown that it is possible to 
put specific constraints on the nature of the interacting modes 
in a couple of variable white dwarfs. 

We have discussed in some detail the relation between the 
frequency locks that can be introduced by resonances, and the 
constancy of the associated modal amplitudes. Conversely we 
have also discussed the amplitude and frequency modulations 
that can appear beyond the immediate vicinity of the resonances. 

In the case of a 2: 1 resonance we have shown on the ba- 
sis of a numerical example the following scenario. When the 
star evolves toward a resonance condition it first undergoes pul- 
sations with rapid, but small amplitude and frequency mod- 
ulations. These modulations increase in size and decrease in 
frequency as the resonance is approached. In the vicinity of 
the resonance the star suddenly jumps into constant frequency 
locked pulsations. 

For rotationally split multiplets the presence of frequency 
lock depends on the ratio of (D/u;) 2 being small (in order of 
magnitude) compared to k/uc (Eqs. 16 and 17). When Q is suf- 
ficiently large or k is sufficiently small one obtains amplitude 
and frequency modulated pulsations. 

The general question of ’mode selection’ also arises when 
there are possible several coexisting but mutually exclusive sta- 
ble pulsational states. Which one does the star choose to pul- 
sate in? In general, as evolution proceeds, the choice of pul- 
sational state is unique. The simplest example is given by the 
nonresonant interaction of two modes ( cf. Buchler & Kovdcs 
1 986). At the instability blue edge the star starts to pulsate in the 
mode which has just become linearly unstable (more precisely 
by mode we mean the nonlinear mode which is an outgrowth 
of the unstable linear mode). It will then stay in this mode until 
this mode itself becomes unstable or bifurcates into another one. 
There can be hysteresis. During the rightward evolution in the 
HR diagram the pulsational state may differ from the one a left- 
ward evolution would give. For completeness we recall that the 
coexisting stable pulsational states which can only be reached 
through a finite perturbation away from a stable pulsational state 
in which the system finds itself. Whether intrinsic turbulent or 
convective fluctuations could be large enough to effect such a 
“hard” bifurcation is unknown. 

The nonlinear asteroseismology that we have discussed here 
is still in its infancy. We have shown how resonances can be 
used to help in the identification of modes. A few interesting 
results can be obtained from a phenomenological confrontation 
between the predictions of AEs and the observations. The major 
quantitative progress that we foresee, though, will come with 
the wedding of a b initio computations of the coupling coef- 
ficients to more comprehensive observations. With respect to 


the latter, great strides have been made with multi-site obser- 
vations (WET, STEPHI, etc.). These collaborative efforts have 
been mostly geared toward linear asteroseismology and it would 
be very useful if they could be more attuned to specifically non- 
linear considerations; in particular concentrate on those w'hite 
dwarfs that have simple spectra, with ostensive resonances, both 
those with apparently constant amplitudes and those w ith mod- 
ulations. 
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